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Abstract. Precision measurements on nucleons provide constraints on the Standard 
Model and can also discern the signatures predicted for particles beyond the Standard 
Model. Knowing the Standard Model inputs to nucleon matrix elements will be 
necessary to constrain the couplings of dark matter candidates such as the neutralino, 
to relate the neutron electric dipole moment to the CP-violating theta parameter, 
or to search for new TeV-scale particles though non-U— A interactions in neutron 
beta decay. However, these matrix elements derive from the properties of quantum 
chromodynamics (QCD) at low energies, where the coupling is strong and thus 
perturbative treatments fail. Using lattice gauge theory, we can nonperturbatively 
calculate the QCD path integral on a supercomputer. 

In this proceeding, I will review a few representative areas in which lattice QCD can 
contribute to understanding the structure inside nucleon and how they can contribute 
to the search for beyond-the-Standard Model physics, with discussions of the difficulties 
and prospects for future development. 
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1. Introduction 

Quantum chromodynamics (QCD) is the theory of the physics that dominates at very 
small scales, the femtoscale (10~ 15 m), where fundamental particles (quarks and gluons) 
are bound inside the building blocks of matter, protons and neutrons. This is far 
smaller than the size of the atoms (10~ 10 m) which form molecules (nanoscale, 10 -9 ) in 
the materials we are familiar with, like a cup of coffee. What do we learn by considering 
physics at this tiny scale? For one, we can understand how fundamental particles make 
up the proton and neutron, and more exotic particles such as hyperons. Further, we 
may try to see how they interact with each other. At the nuclear level, we can learn 
how the residual strong force binds nucleons together into the nuclei that form the 
cores of all atoms, that power the stars and perhaps one day fusion power plants. An 
understanding of how nuclear structure arises from fundamental physics may hold the 
key to the question of fine tuning (a mystery relevant to any carbon-based life forms). 
The interactions between hadrons at small scale also impact objects at the astrophysical 
scale, such as neutron stars. The strengths of these interactions can be input into 
descriptions of nuclear matter in simulations of the conditions in the center of neutron 
stars, where incredible pressure and density could allow exotic forms of matter to exist. 

Quantum chromodynamics is the theory of the color force. It describes the strong 
interactions between quarks and gluons using an SU(3) gauge theory. Given a QCD 
action S and interesting observables whose properties are described by an operator O, 
one can compute the physical quantities of interest using a path integral (integrating 
over all possible configurations of gluonic and fermionic fields A, ip, ■0): 



One of the interesting properties of quantum chromodynamics is confinement. That is, 
we never see free quarks in nature; rather, all we see are composite particles containing 
quarks, called hadrons, such as the lightest hadron, the pion, and the proton and neutron 
of ordinary nuclei. At large energy, the QCD coupling is small. We can simply make an 
expansion in terms of the coupling, and it converges well. In fact, the calculations can 
be easy enough (at low order) to give to your graduate field-theory class as homework. 
However, at low energy the strong coupling becomes larger and larger, and, by say 
the scale of neutron beta decay, perturbation theory becomes poorly convergent. At 
low energy, the theory becomes nonperturbative and even just the vacuum of QCD is 
incredibly complicated. The QCD vacuum teems with topological charge and the results 
of the broken quark chiral symmetry. Figure [T] shows a few timeslices of the spatial 
distribution of topological charge in an example QCD vacuum configuration. Unlike 
the classical vacuum, just to describe this is incredibly complicated and a task that is 
nearly impossible analytically. Therefore, nonperturbative approaches are essential to 
exploring QCD physics at this energy scale. In order to study the interesting physics in 
the low-energy regions, in the 1970s Kenneth Wilson proposed to discretize space and 
time in the path integral and to work in Euclidean space. Thus, lattice QCD (LQCD) 
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Figure 1. A selection of timeslices showing the topological charge of the vacuum 
calculated on one gauge-field configuration on a lattice with spacing approximately 
0.1227 x 10 -15 meters with 16 lattice sites in each spatial direction. The evolution 
of the vacuum shown in this image occurs over 5.8 x 10~ 25 seconds. The red regions 
correspond to instantons with positive topological winding number, while the blue 
correspond to instantons with negative winding number. 

was born. It contains two scales that are absent in continuum QCD, one ultraviolet (the 
lattice spacing a) and one infrared (the spatial extent of the box L). The problem now 
only involves a finite number of degrees of freedom and can be put on a computer for 
numerical integration. 

Of course, if computational resources were not an issue, we could solve all the QCD 
problems by brute force as long as we could write down the path integrals. However, it 
is not possible within a human lifetime to do so, so we have to examine the problem and 
solve it more cleverly. The number of degrees of freedom for a fermion here is the spatial 
volume times the number (3) of colors, times the number (4) of spins. The number of 
spatial points varies, for example, a lattice volume with 64 3 x 96 has rank 300M, and it 
cannot be fit onto a few workstations. A computational cluster with thousands of cores 
is necessary. Lattice practitioners use systems ranging from small-size clusters hosted at 
universities to national supercomputer centers, such as Kraken at NICS, or CPU-GPU 
hybrid systems, such as Edge at LLNL. 

In most cases, we spend the majority of computational hours inverting the Dirac 
operator. Due to the breaking of continuous rotation symmetry, multiple versions of 
the Dirac operator can be written on the lattice that describe the same physics in the 
continuum. Such a discrete Dirac operator is a sparse and structured matrix, properties 
that we can take advantage of to develop smarter solvers for our problem. To illustrate 
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Figure 2. (left) Two-dimensional illustration of the spin structure of the Wilson 
Dirac operator. Gray is zero; blue is real and positive; red is negative; magenta and 
green are positive and negative imaginary, (right) Spectral plot for a two-dimensional 
Wilson Dirac operator. The blue circles indicate the results from free-field, while the 
olive dots have U(l) interacting gauge fields. The smaller the m q , the worse the matrix 
condition number. 

structure of our operator, we show a 2-dimensional version with U(l) gauge fields (in 
reality, we use a much more complicated 4D, SU(3) theory); this particular form is the 
Wilson Dirac operator: 

^ son = (m+d) <W-£ \ K 1 + 7„)^<W + (1 - 7^,A*'+a] -(2) 

The left-hand side of Fig. [2] shows the spin structure of the two-dimensional version of 
Eq. [2| We can see that it is highly structured along diagonals, with non-zero terms in 
regular places, but mostly zero. The eigenvalues of this operator in the complex plane are 
shown on the right-hand side of Fig. [2] When interactions are turned off, the eigenvalues 
are these blue circles, while the olive points are the eigenvalues with thermalized U(l) 
gauge fields included. As the quark mass is decreased, entire eigenspectrum is shifted 
to the left; thus, the lowest eigenvalues approach zero and the matrix condition number 
diverges. This means that our inversion algorithms will run very slowly. Only recently 
have ensembles with physical pion masses become available; most of the time, we 
do the calculation with multiple heavier masses, and a procedure for extrapolation is 
necessary. Taking advantage of the structure of the operator, many techniques have been 
developed to solve the Dirac equation faster: eigenvector deflation, multi-mass solvers, 
multigrid acceleration, stochastic sources. Some such techniques have been examining 
the memory /communication problem, which will be needed to run QCD effectively on 
systems including GPUs PQ. 

Once we have D~ l , we can use it to calculate Green functions from which we will 
extract physical quantities. For example, we can calculate the pion mass by creating a 
pion at a certain point in space and time and annihilating it at another. In practice, we 
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"contract" the inverse Dirac matrix with some spin/color matrix: 

CV(t) = <7r(0)|7r(t)> = (u l5 d(0) d l5 u(t)) = (D-\0,t)(D-\0,t)y) . (3) 

After momentum projection to the rest frame, the resulting Green function declines 
exponentially in the time direction with exponent equal to the mass. 



2. Probing inside Nucleons 

To "see" inside a hadron, we need a stronger "microscope" than the typical optical ones; 
to probe inside the proton requires, for example, a GeV-scale electron beam elastically 
scattering with a nucleon. This creates interaction with the nucleons through exchange 
of photons; in QCD, it can be simulated as a vector current to probe electromagnetic 
properties. We still create the particle and annihilate it at another coordinate of 
spacetime; in addition, we now insert a vector current to probe the nucleon's structure 
in the following general form: 

= (N(x 2 )J(x 1 )N(x )), (4) 

where 

J(x) = ^2e f q f (x)0(x,y)q f (y), (5) 
/ 

and we sum over all flavors / G {u, d}; O can be either simple gamma matrices or their 
products with some covariant displacement function. O = j/j,S Xiy is used for simulating 
the elastic scattering process. Such a current can connect directly to the quark fields in 
the baryon operator or it can interact with vacuum quark loops. The available momenta 
in the system are limited by the lattice's boundary conditions, the box size and lattice 
spacing. By varying the inserted operator J, we can study various properties of nucleons. 
In this section, we will present a few examples on nucleon axial charge, proton spin and 
neutron transverse density. 



2.1. Nucleon Axial Charge qa 



The most-calculated nucleon matrix element is the nucleon axial charge, Qa- It is one 
of the most fundamental properties of the nucleon; experimentally, it can be measured 
though neutron beta decay. It is also an important parameter that determines the 
rate of proton-proton fusion, the dominant hydrogen-burning process inside the sun. It 
can be extracted from the three-point Green function, as shown earlier, but with an 
axial-vector current inserted: 

G P (q 2 ) 



(N\A„\N)=m(p') 



2M, 



A' 



un(p), 



(6) 



where Ujy are nucleon spinors and q = p' — p is the momentum transfer. The nucleon 
axial charge is defined as qa = G A ^ d (q 2 = 0). The left-hand side of Fig. [3] shows a 
collection of such calculations as a function of quark mass (parametrized by pion mass), 
from all Nj = 2 + 1 (degenerate up/down and a strange flavors in the sea) lattice 
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Figure 3. (left) Summary of the lattice Nf = 2 + 1 results on gA as functions of 
[2j [3l HI [3 [6] . (right) The world average of gA using Nf = 2 and 2 + 1 dynamical 
lattices results [3l [TJ El HI O El [10l [11]. 



calculations. Most of the results here only display statistical errors, but nonetheless, we 
can see that the various versions of the Dirac operator generate similar results for this 
quantity. Taking an average over all the dynamical results extrapolated to the physical 
pion mass, we obtain an LQCD g& of 1.16(3) with about 3% error, as shown in the 
right-hand side of Fig. |3j 

2.2. Proton Spin 

Another fundamental question is how the proton's spin 1/2 is distributed among its 
constituents. The most naive intuition is that the three quarks carry the spin. Their 
contributions can come from their intrinsic spin and orbital angular momentum. We use 
the set of covariant operators suggested by Ji [12] and calculate. Again, we show results 
as a function of light-quark mass (pion mass). Figure [4] shows a selection of example 
results calculated by different collaborations (Nf = 2 QCDSF [13] and Nf = 2 + 1 
LHPC [I]) who use different Dirac operator interpretations (clover and domain- wall 
fermions on staggered sea, respectively). These results are renormalized in the MS 
scheme at a scale of 2 GeV. The extracted values are the stars on the left. As it turns 
out, the quarks contribute less than 50% of the total nucleon spin; the majority of it 
comes from gluons, which is quite a surprise. However, both of the calculations have 
ignored the "disconnected" diagram; that is, the inserted operator constructed such that 
its quarks only connect with the valence nucleon quarks through gluons. Calculating 
such contributions requires more computational resources to get clean signal, but they 
could contribute up to 20-30% of the total quark contributions [H] , A few groups have 
made remarkable progress on this type of contribution [151 EH El H]> as wen as 011 
direct gluon contributions to the spin [19|. 



Lattice QCD for Precision Nucleon Matrix Elements 



7 



0.4 
0.3 
0.2 
0.1 


-0.1 



-, — I — , — I — , — I — , — I — , — I — , — I- 
LHPC, QCDSF 

* * * * 



LHPC, QCDSF 

i U ^ — , 



KI!' +d l2 



r u+d 



o 



0.2 



0.4 0.6 

M 2 (GeV 2 ) 



0.8 



0.5 
0.3 
0.1 
-0.1 
-0.3 





i . i . i 


* 


i 








AX"/2 : 




▼ ▼ 






* <• 


AI. d /2 
i ... i 


\- 


i 







0.2 0.4 

M 2 (GeV 2 ) 



0.6 



Figure 4. (left) Total quark contribution to proton spin and and orbital angular 
momentum as a function of pion masses by QCDSF and LHPC. The total orbital 
angular momentum from quarks is about 6%. (right) Decomposition of the individual 
quark contribution to spin and angular momentum by LHPC. The individual up 
and down quark contribution to angular momentum is about 36% in magnitude but 
opposite in sign; as a result, the total quark contribution is mostly canceled, leaving 
total around 6%. The up-quark spin contribution to the total proton spin is about 
82% but the down-quark contributes about half of the size in the opposite sign; leaving 
total quark spin contribution to proton spin of 40% or so. 



2.3. Transverse Densities 

Since LQCD gives us access to a complete sample of the QCD vacuum in spacetime, we 
can also investigate the distribution of charge density within the nucleon. To avoid 
relativistic ambiguity due to the transfer of momentum by electromagnetic probes, 
we should examine the transverse charge density (in directions perpendicular to some 
impact plane, as an ultrarelativistic electron would see) for a polarized nucleon. Firstly, 
we need to obtain the nucleon vector-current matrix element over a large range of 
momenta and from these obtain the electromagnetic form factors {Fi^{Q 2 )) ■ Then by 
taking a sort of Fourier transform, we can see how the transverse charge density is 
distributed as a function of the impact distance b in a polarized nucleon [201 [2D 122] : 

PTib) = r^-J,(bQ)F,{Q 2 ) + sin(0) rf£l Ji(bQ)F 2 (Q 2 ), (7) 
Jo ^ Jo lnM N 

where Jo,i are Bessel functions. We can perform this integral numerically, using the 

lattice Fi 2 {Q 2 ) obtained by extrapolating our fit form to the physical pion mass. Taking 

the large-Q 2 form-factor data from Ref . [23] , Figure [5] shows the results for the neutron 

in one dimension (right, shown as the blue band here) and the two-dimensional impact 

plane (left) using our lattice inputs. There are positive and negative charges surrounding 

the center, which in the neutron sum to zero. The dashed line on the left-hand side 

is a similar application of Eq. [7] using experimental neutron form factors, which are 

only known up to 1.5 GeV 2 ; thus, a majority of the form-factor inputs are based on 

extrapolation to the larger-Q 2 region. Note that the asymmetry in the distribution for 
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Figure 5. The transverse charge densities in a polarized neutron with b XiV (left) 
and b (right) in fm. In the left-hand figure, black indicates near-zero values, and 
purple, orange and white are increasingly positive. In the right-hand figure, the blue 
band indicates the densities from 2+1-flavor lattice calculation [53], while the dashed 
line is a parametrization interpolating and extrapolating from available experimental 
data [25] . (Note that the neutron experimental data is extrapolated to larger-Q 2 
region). 

a polarized nucleon is due to the relativistic effect of boosting the magnetic moment of 
the baryon. This induces an electric dipole moment that shifts the charge distribution. 

3. Applications beyond the Standard Model 

Can we learn about physics beyond the Standard Model (BSM) with LQCD nucleon 
matrix elements? There are many opportunities for LQCD to provide important inputs 
to probe BSM physics with nucleons. For example, the presence of new particles beyond 
the Standard Model at the TeV scale can be probed through detection of new scalar and 
tensor interactions in neutron beta decay. New ultra-cold neutron (UCN) experiments, 
such as the UCNb/B experiment at Los Alamos National Lab (LANL) or one proposed 
to run at Oak Ridge National Lab (ORNL) in 2015 can provide high-precision low-energy 
data to constrain the allowed range of new-particle masses. The strange contributions 
to the proton's scalar and spin densities are important inputs for spin-independent and 
-dependent dark-matter cross-section measurements. 

3.1. Non-V—A Interactions in Nucleon Beta Decays 

Firstly, we can try to look for non-Standard Model contributions in precision neutron 
(nuclear) beta-decay measurements. The neutron beta-decay Hamiltonian contains the 
Standard-Model V— A current for the leptons and the quarks, and it may also contain 
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terms corresponding to new BSM physics: 



tfeff = G F Jj* X J<™* + ^efSM^lept x ^ 



quark 



) 



(8) 



where Gf is the Fermi constant, Jv-a is the left-handed current of the indicated particle, 
and the sum includes operators with novel chiral structure. So in the context of our 
theory, new operators will enter with the coefficients e that are related to the TeV scale 
of the particles. The leptonic part is understandable using analytic techniques, but the 
quark operator in the context of the nucleon will introduce some unknown coupling 
constants: 



which are nonperturbative functions of the nucleon structure, described in the SM by 
QCD. Any deviation from the SM V— A current coming from new scalar and tensor 
interactions in the effective theory will require knowledge of the couplings gs and qt to 
understand. For more details about experimental and theoretical work on this subject, 
we refer readers to Ref. |26j, and references within. 

The tensor charge g? has been studied a few times in the past using Nf = 2 + 1 
dynamical ensembles in lattice QCD, as summarized in the left-hand side of Fig. [6j 
We take the 2+1-flavor results by the RBC [27] and LHPC [28], and recently by 
PNDME [2S], and make a global plot as a function of pion mass squared, gs, on 
the other hand, has not been much studied; a few collected calculations [29] is shown 
on the right-hand side of Fig. [6] using various Nf = 2 + 1 lattice actions. All of these 
results, by different collaborations, using different fermion actions, (but at the same 
lattice spacing) are generally in good agreement, although the errorbars shown here 
contain only statistical error. We further globally analyze all the available lattice data 
(with a upper pion- mass cut at 550 MeV), including (excluding) the data from PNDME 
(whose lightest pion mass is around 220 MeV), shown as the purple (orange) band in 
Fig. [6j We use the chiral formulation given in Ref. [30] and a linear ansatz for tensor 
and scalar charges, respectively, to extrapolate to the physical pion mass. We see that 
the PNDME points greatly constrain the uncertainty due to chiral extrapolation in both 
cases and obtain g^ QCD = 0.95(5) and g^ CI) = 0.69(9). The LQCD values are better 
determined than other theoretical estimations (from different model approximations), 
which give rather loose bounds on these quantities; for example, gs is estimated to be 
between 0.25 and 1. 

Finally, these low-energy matrix elements, the tensor and scalar charges, can be 
combined with experimental data to determine the allowed region of parameter space 
for scalar and tensor BSM couplings (denoted e). Using the gs,T from the model 
estimations and combining with the existing nuclear experimental data[§J we get the 

§ For example, nuclear beta decay + — > + transitions and other processes, such as /? asymmetry 
in Gamow- Teller 60 Co, longitudinal polarization ratio between Fermi and Gamow- Teller transitions 
in 114 In, positron polarization in polarized 107 In and beta-neutrino correlation parameters in nuclear 
transitions. 



9t = (n\ua^ u d\p), g s = (n\ud\p), 



(9) 
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Figure 6. Global analysis of all Nf = 2 + 1 lattice calculations of gr (left) and 
gs (right). The leftmost points are the extrapolated values at the physical pion mass. 
The two bands show different extrapolations in the pion mass: with and without the 
lighter PNDME collaboration points from Ref. [25]; the final uncertainty due to the 
chiral extrapolation is significantly improved. 

constraints shown as the outermost band on the left-hand side of Fig. [7j An ongoing 
UCN experiment is studying neutron beta decay at LANL to look for deviations from 
the SM in the Fierz term and the neutrino asymmetry parameter to the level of 10 -3 by 
2013. Combining those expected data and existing measurements, and again, using the 
model inputs of gs,T, we see the uncertainties in es,r are significantly improved. Finally, 
using our present lattice-QCD values of the scalar and tensor charges, combined with 
the expected 2013 precision of experimental bounds on deviation of these neutron-decay 
parameters from their SM values, we found the constraints on e$,r are further improved, 
shown as the innermost region. These upper bounds on the effective couplings es,T 
correspond to lower bounds for the scales A$,t at 2.9 and 5.6 TeV, respectively, for new 
physics in these channels. 

How do the constraints from high-energy experiments compare? As demonstrated 
in Ref. |26J, neither CDF nor DO is sufficient to provide useful constraints. We can 
estimate the es,r constraints from LHC current bounds and near-term expectations 
through effective Lagrangian 

C = -^V ud (ud)(eP L u e ) - ^V ud (ua^d)(ea^P L u e ) } (10) 

Jig Arp 

where t)s,t — ±1- By looking at events with high transverse mass from CMS/ATLAS in 
the eu + X channel and comparing with the SM W background, we estimated 90%-C.L. 
constraints on r]s } T based on the current LHC run, y/s = 7 TeV C = 1 fb" 1 (the green 
line) and for a near- future run y/s = 7 TeV £ = 10 ftT 1 (the purple dashed line) on the 
left-hand side of Fig. [7j More details can be found in Ref. [26] . 

3.2. Dark Matter Cross-Sections 

Recent evidence suggests that the composition of dark matter must be cold, with 
velocities that are non-relativistic, giving it a Compton wavelength above the width 
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Figure 7. (left) es-er allowed parameter region using different experimental and 
theoretical inputs. The green and purple dashed lines are the constraints from LHC 
current bounds and near-term expectation, (right) Summary of dynamical LQCD 
calculations of the nucleon sigma term m s (ss) in MeV. The (pink) band indicates the 
global average after weighting individual points according to their systematics. 



of a proton. One promising candidate for dark matter is the supersymmetric particle 
called the neutralino. If dark matter is neutralinos, they will interact with nuclei 
at low-energy. A typical value for the cross-section in the Constrained Minimal 
Supersymmetric Standard Model (CMSSM) is about 1 zeptobarn (10 -49 m 2 ). The 
strange-quark contribution to the proton scalar density and spin are important inputs 
for spin-independent and -dependent cross-sections. In this proceeding, we will focus on 
the LQCD calculation of a s = m s (N\ss\N). 

Many groups have developed techniques to investigate a s . One can directly 
calculate the matrix element which involves strange-quark loops and proton correlators. 
Unfortunately, such a direct calculation is noisy, so various techniques have been 
developed to improve the signal (e.g. JLQCD, QCDSF, BU, Engelhardt) [HH E2J LEI EH1 
EH [35], |36] . Alternatively, one could use the Feynman-Hellman Theorem by studying how 
the nucleon mass varies with the strange-quark mass, taking the numerical derivative 
dM^jdm s either by direct SU(3) fitting to the baryon masses (e.g. Young/Thomas) or 
by reweighting the strange part of action (e.g. Jung, MILC) [371 EH EH] . The values of a s 
from various calculations are shown on the left-hand side of Fig. [7} We perform a global 
fit to current dynamical lattice data (weighted by lattice spacing, lightest m n , dynamical 
strange and other quality factors) and obtain <jLQcp_ 43(g) MeV. The value can be input 
into theoretical calculations of cross-section constraints from various models; an example 
of the application to the CMSSM can be found in Ref. [30] . 
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4. Road to High Precision 

Nucleon matrix elements are generally more computationally demanding than mesonic 
matrix elements, since the statistical noise grows with Euclidean time t as e ( M N-3M 7T /2)t 
for each nucleon in the system. Such characteristic signal-to-noise behavior has 
been studied and well-demonstrated up to small nuclei using high statistics by 
NPLQCD [H], H2J H3j H3]. Take g^, for example; the present uncertainty (including 
the systematics) is around 6-10% pEJ EJ EJ HDj- If we are to target statistical error 
for qa at a level of 1% or smaller, we need at least an order-of-magnitude boost in 
computational resources. 

Anther obstacle to obtaining high-precision nucleon matrix elements is to have 
control over the systematic uncertainties of the LQCD calculation. The lattice scales, 
the lattice spacing and box size, must be taken a — > and L — > oo, or we must at least 
show that the quantities of interest are independent of these unphysical scales. The 
nearest excited state of the nucleon, the Roper resonance, is only few hundred MeV 
away, unlike the n', which is almost an order of magnitude heavier than the pion. As a 
result, it can potentially create significant contamination to the wanted nucleon matrix 
element. One way to control this systematic is to increase the source-sink separation in 
the matrix-element calculation, since the excited-state contamination is exponentially 
suppressed; however, this will also decrease the signal-to-noise ratio of the desired matrix 
element as the separation enlarges. Refs. [21 [23] have explicitly taken account of the 
excited-state contribution to the nucleon matrix element in their analysis. 

Chiral extrapolation of nucleon quantities can also introduce large systematics, 
especially when the formulation contains significant contributions from chiral 
logarithms. Figure [6] clearly demonstrates how an extrapolation could be improved 
by the inclusion of calculations at lighter pion mass. Such extrapolations also raise the 
questions of the validity of chiral perturbation theory over the range of pion masses 
used, how a particular expansion scale adopted converged, and whether one should use 
expansions with SU(3) or SU(2) flavor symmetry. In all cases, having light pion masses 
is essential and many collaborations already have begun generating gauge ensembles at 
the physical pion mass, which will soon remove extrapolation systematics. 

To compare lattice results with continuum ones, we also need to renormalize 
the lattice operators used in Eq. [4] and match to the continuum MS scheme 
(running to 2 GeV scale). Unlike the continuum case, calculating the lattice-operator 
renormalization factor in perturbation theory is complicated and slowly converging. To 
reduce the systematic associated with renormalization, it is commonly agreed that using 
a nonperturbative approach (such as Schrodinger Functional or RI/MOM scheme) on 
the lattice should work best. To extract high-precision nucleon matrix elements, we 
need to improve the accuracy of nonperturbative renormalization as well as working out 
the continuum matching to several-loop order for each operator of interest. 

With the ever increasing computational resources made available to us, we can 
overcome the challenges that occur in nucleon systems to improve the uncertainties on 
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nucleon matrix elements. Furthermore, extrapolation to the physical light-quark masses 
is likely to be brought under control in the next few years, as ensembles of lattices begin 
to be generated with physical u, d (and s and c) quark masses, which should greatly 
reduce the systematic uncertainties. Other systematics, such as finite-volume effects, 
renormalization and excited-state contamination can be reduced by improved algorithms 
and by increasing the computational resources devoted to the calculations. 

To sum up, the name of the game is precision, and lattice QCD is now entering an 
exciting era where we can provide precision input directly from the first principles of 
the Standard Model. This is partially thanks to the global increase in computational 
resources available to us, and also due to improved algorithms, which have greatly 
improved our abilities to calculate quantities that were not accessible before. There are 
many quantities that we can contribute using supercomputers. We see various different 
lattice actions and groups with independent calculations provide consistency checks. 
We have reproduced well measured experimental values and provided better SM values. 
We can make predictions for quantities that have not or could not be measured by 
experiment. Lattice QCD is essential in cases where experiments are limited, such as 
for the scalar charge or hyperon axial couplings. 

New experiments at both the high-energy frontier and the precision low-energy 
frontier are now enabling us to probe BSM physics. But to do so, we need precision 
quantities from Standard Model, such as gs,r or the strangeness in proton, to constrain 
various possible BSM models. By continuing to improve the precision of our nucleon 
matrix elements, we can play a considerable role in this search. We will also be able 
to identify high-impact contributions that lattice techniques can make. These could 
be for quantities which are difficult or impossible to measure in experiment and where 
models introduce unacceptable levels of uncertainty. By determining what needs to be 
computed and delivering precise values from first principles, lattice QCD will solidify 
its role in particle and nuclear physics. 
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